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Introduction 


Abstract 


The Interior Highlands, in southern North America, possesses a distinct fauna with nu- 
merous endemic species. Many freshwater taxa from this area exhibit genetic structuring 
consistent with biogeography, but this notion has not been explored in freshwater snails. 
Using mitochondrial 16S DNA sequences and ISSRs, we aimed to examine genetic struc- 
turing in the Pyramid Elimia, E/imia potosiensis, at various geographic scales. On a broad 
scale, maximum likelihood and network analyses of 16S data revealed a high diversity of 
mitotypes lacking biogeographic patterns across the range of E. potosiensis. On smaller 
geographic scales, ISSRs revealed significant population structure, even over the distance 
of a few hundred meters. Unlike other freshwater mollusks like mussels, FE. potosiensis 
showed no evolutionary patterns relating to biogeography. The species does show popula- 
tion-level genetic structure, which may have implications in conservation efforts. 


riods of glaciation and sea level rise. Analyses of Inte- 
rior Highland fauna demonstrate a specific relationship 


The Interior Highlands, separated by the Arkansas River 
Valley into the Ozark and Ouachita regions, is a major 
geographical feature of southern North America and in- 
cludes the Ozark and Ouachita Mountains (Fenneman 
1928, Mayden 1985). As a result of the geological pro- 
cesses that formed it, the Interior Highlands has a distinct 
fauna with high instances of endemism and drainage pat- 
terns that have dictated diversity (Mayden 1985, Mayden 
1988, Matthews and Robison 1998, Austin et al. 2004, 
Bonett and Chippindale 2004). The Interior Highlands 
also served as a refugium for many species during pe- 


between southern Ozark fauna (from Arkansas River 
tributaries) and Ouachita Mountain fauna (Mayden 1988, 
Crandall and Templeton 1999, Berendzen et al. 2003). 
For some aquatic species, the Arkansas River may have 
formed a barrier for those adapted to cold, fast moving 
waters (Turner et al. 1996, Bonett and Chippindale 2004), 
thus isolating populations. Genetic structuring at narrow 
and broad geographic scales have been observed in fresh- 
water bivalves from the Interior Highlands (Inoue et al. 
2013, Chong et al. 2016) but has not been explored in 
freshwater snails in the same region. 
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With nearly 170 recognized species (Johnson et al. 
2013), the family Pleuroceridae is the second largest 
group of freshwater snails in North America behind 
only the Hydrobiidae. As part of the 700 snail species in 
North America, pleurocerids are increasingly imperiled 
by river regulation, habitat loss, poor water quality, re- 
duced water quantity, and invasive species (Johnson et 
al. 2013). Only 22% of all freshwater snail species are 
stable in a conservation context, and recent extinctions 
support these rankings (Sada and Vinyard 2002, Hershler 
et al. 2007, Johnson et al. 2013). Thirty-three pleurocerid 
Species are extinct, and of all federally endangered and 
threatened species, 55% of the listed snails are fresh- 
water species including eight pleurocerids (Johnson et 
al. 2013). Only six recognized species occur west of the 
Mississippi River (NatureServe 2017), with three of them 
occurring in the Interior Highlands. The Pyramid Elimia, 
Elimia potosiensis (Lea, 1841), is a morphologically 
variable species found in springs and rivers throughout 
the Ozark and Ouachita highlands of Arkansas and Mis- 
sourl, its range stretching westward into eastern Kansas 
and Oklahoma (Goodrich 1939, Tiemann and Cummings 
2007). Unlike many other pleurocerids, populations of 
E. potosiensis are often large and connected to one an- 
other within stream flows (Gordon 1980, Gordon 1982), 
and the species is secure across its range (global heritage 
rank of G5; NatureServe 2017). Given that freshwater 
taxa, including mollusks, in the Interior highlands exhib- 
it genetic structure, and that many freshwater snails are 
of conservation concern, we aimed to answer three re- 
search questions. First, on a broad geographic scale, does 
E. potosiensis represent a single clade showing genetic 
structuring by population and river drainage? Second, on 
narrower scales, do populations show any degree of ge- 
netic structure within single streams of closely connected 
waterways? Finally, what impact can our findings have 
on pleurocerid conservation? 

We utilized two forms of genetic data for our study, 
mitochondrial DNA sequences at the species level and 
inter-simple sequence repeats (ISSRs) at the population 
level. The bulk of freshwater gastropod genetic studies 
employ one of two mitochondrial DNA fragments de- 
rived from either the 16S or cytochrome oxidase c sub- 
unit I gene. While these sequences may show utility in 
population and species-level studies, results are mixed 
in freshwater cerithioideans including pleurocerids (e.g. 
Kohler and Deein 2010, Miura et al. 2013, Kohler 2016). 
Widely divergent mitochondrial haplotypes can be found 
in single populations and species that can cloud questions 
of monophyly and relatedness (Dillon and Frankis 2004, 
Whelan and Strong 2016), and no consistent explanation 
has been offered as to why these divergent haplotypes ex- 
ist or persist (Whelan and Strong 2016). Not all studies, 
however, report this issue (Lydeard et al. 1998, Minton 
and Lydeard 2003, Minton and Savarese 2005). ISSRs are 
fragments of DNA separating neighboring microsatellites 
used as genetic markers amplified using polymerase chain 
reaction (PCR) (Zietkiewicz et al. 1994). Primers for 
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ISSR amplification are complementary to microsatellite 
sites, bind to them using a one to three nucleotide anchor 
sequence on the 5’ or 3’ end (Bussell et al. 2005), and am- 
plify the ISSR in between (Culley and Wolfe 2001). This 
approach does not require prior genome sequence infor- 
mation, and leads to banding patterns containing multi- 
ple loci and high levels of polymorphism (Tsumura et al. 
1996). ISSRs have provided useful population data in a 
variety of organisms including plants (Lisek and Rozpara 
2010), insects (Vijayan et al. 2006), arachnids (Machk- 
our-M’ Rabet et al. 2009), and birds (Haig et al. 2003), but 
have not been utilized extensively in gastropods (Dong 
et al. 2011, Snegin 2014). 


Materials and methods 


Broad scale phylogenetic analyses 


We collected live E. potosiensis from six river drainag- 
es throughout the Ozark and Ouachita highlands (Fig- 
ure 1). We followed a standard CTAB/chloroform protocol 
(Saghai-Maroof et al. 1984) for DNA extraction. We ampli- 
fied a 500 bp portion of the 16S ribosomal subunit (primers 
l6sar and 16sbr [Palumbi, 1996]) using Qiagen Taq PCR 
Core Kit and the following thermal cycling conditions: 
an initial denaturation cycle of 95°C for three minutes; 40 
cycles of 95°C for 35s, 44°C for 45s, 72°C for 45s; and a 
final five-minute extension period at 72°C. We sequenced 
gel-purified products on a Beckman CEQ8000 automated 
sequencer using Sanger dideoxy sequencing. We used Po- 
part (Leigh and Bryant 2015) to calculate nucleotide diversi- 
ty and Tajima’s D (Tajima 1989), and to visualize haplotype 
diversity in E. potosiensis by generating a TCS (Clement et 
al. 2002) network. We then combined our E. potosiensis se- 
quences with other pleurocerid 16S data taken from NCBI 
(Table 1) and aligned them using the Muscle algorithm 
(Edgar 2004) with default settings. We used Gblocks 0.91b 
(Castresana 2000) with default settings to remove poorly 
aligned positions from the data, and used the ModelFinder 
(Kalyaanamoorthy et al. 2017) function of Iq-tree (Nguyen 
et al. 2015) to select the appropriate substitution model by 
Bayesian information criterion. We then analyzed the data 
under maximum likelihood in Iq-tree and estimated branch 
support using 10,000 ultra-fast bootstrap replicates (Minh et 
al. 2013). We also tested the hypothesis that E. potosiensis 
sequences from each separate river drainages represented 
separate monophyletic groups. We used the approximately 
unbiased test (Shimodaira 2002) with 1000 RELL bootstrap 
replicates implemented in Iq-tree to test whether topologies 
were significantly different (p<0.05). 


Narrow scale population genetic analyses 


We employed ISSRs as population markers on two differ- 
ent spatial scales. In our first study, herein referred to as 
the ‘small’ study, ten individuals from each of 12 E. po- 
tosiensis populations were collected (Figure 2, white cir- 
cles). Four populations were in a spring run flowing into 
Walnut Creek in the Arkansas portion of the Ouachita 
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Figure 1. Map of E. potosiensis collection localities, color coded by river drainage, used in the phylogenetic analysis. Detailed 


locality information is in Table 1. 


National Forest. Populations were sampled at the mouth 
of the spring run (S4) and in 50 m intervals upstream (S3- 
S1); at the mouth, the spring run flowed over a 1 m rock 
ledge into Walnut Creek. Populations from Walnut Creek 
were sampled as well, using the confluence point with 
the spring drainage (C3) as a starting point. Populations 
were sampled in 50 m intervals upstream (C2 and Cl) 
and downstream (B1—B5, below where the creek flows 
under the road through a culvert). 

We followed the ‘small’ study with a second ‘large’ 
study. We collected 50 E. potosiensis from each of four lo- 
calities (Figure 2, gray circles). Three sites were located in 
the Walnut Creek drainage in the Arkansas portion of the 
Ouachita National Forest. The first (ARK1) was a small 
spring feeding the main channel of Walnut Creek. The sec- 
ond and third were in Walnut Creek, 325 m upstream and 
350 m downstream of where the ARK1 spring run entered 
the creek. The final site was Sallisaw Creek, Sequoyah 
County, Oklahoma near Marble City. For both studies, sam- 
ples were stored in 95% ethanol or frozen at -20°C. 


We isolated genomic DNA as before, purified it on 
silica filters (UltraClean PCR Clean-up Kit, MoBio), and 
diluted it with sterile water to 50 ng/ul concentration. We 
amplified ISSRs by PCR in 50 wl volumes with the Go- 
Taq PCR Core System I reagents (Promega) at the man- 
ufacturer’s recommended concentrations. The cycling 
profile consisted of an initial denaturation cycle of 95°C 
for three minutes, 30 cycles of 95°C for 30s, annealing 
for 30s, and 72°C for 60s, followed by a ten-minute ex- 
tension period at 72°C. We used one primer in each ISSR 
reaction (Table 2) and optimized each annealing tem- 
perature. For the ‘small’ study we visualized ISSRs by 
running 10 ul of PCR product on 2% agarose gels in TBE 
at 100 V for 90 minutes. For the ‘large’ study, 20 ul of 
PCR product was run on 10% polyacrylamide-TBE gels 
at 110V for two hours. In both, a 0.1—1 kb size ladder was 
loaded, and we stained gels 1% ethidium bromide then 
de-stained in deionized water. Gels were photographed 
(Bio-Rad Chemi XRS) and processed (Bio-Rad Quantity 
One) before analysis. We coded individual bands if they 
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Table 1. Locality and NCBI accession numbers for sequences used in this study. 


Taxon Accession | Locality Identifier | Reference 
Elimia alabamensis U73761 Lydeard et al. 1997 
E. caelatura AF100988 Holzangel and Lydeard 2000 
E. catenaria FJ471493 Strong and Kohler 2009 
E. clenchi FJ471492 Strong and Kohler 2009 
E. comalensis KU052563 _| Salado Creek at Interstate 35, Salado, Bell County, TX | | new 
E. crenatella U73762 ee Lydeard et al. 1997 
E. cylindracea U73765 SS Lydeard et al. 1997 
E. doolyensis DQ311118 Lee et al. 2006 
E. haysiana U73763 Lydeard et al. 1997 
E. hydei U73764 Lydeard et al. 1997 
E. interrupta AY010521 Lydeard et al. 2002 
E. laqueata KU052565__| Green River, KY new 
E. livescens | DQ311116 | [Lee et al. 2006 
E. livescens 2 KU052564__| French Creek, PA | | new 
E. melanoides AF540003 | | Minton et al. 2003 
E. olivula U73766 Lydeard et al. 1997 
KT988910 7? 
Alum Fork Saline River, AR 34.67310N, 92.79920W 2 new 
KT988911 
KT988965 _| Illinois River, AR 36.10320N, 94.34500W 14 new 
KT988932 
LK War Eagle Creek, AR 36.12100N, 93.69340W 45 new 
KT988962 
KT988940__| Otter Creek, AR 36.22380N, 92.25190W new 
KT988923 
KT988950__| Blanchard Springs, AR 35.95680N, 92.13960W ne new 
KT988951 
KT988956 ae 
Spring River, AR 36.31530N, 91.49080W 123 new 
KT988957 
KT988964 _| Mill Creek Spring, AR 36.05720N, 91.60890W new 
KT988916 
KT988922 
KT988943__ | Mammoth Spring, AR 36.49580N, 91.53320W 154-1 | new 
KT988944 
KT988945 
KT988933 
KT988934 | Warm Fork Spring, AR 36.49580N, 91.53320W 154-2 | new 
KT988929 
KT988954 
E. potosiensis —————— Mulberry Fork Little Red River, AR 35.74210N, 92.33380W 165 new 
KT988955 
KT988939__ | Mulberry River, AR 35.62275N, 93.91023W 186 new 
KT988958 
Cossatot River, AR 34.37950N, 94.23680W 221 new 
KT988967 
KT988963 _ | Tributary to South Fork Ouachita River, AR 34.51870N, 93.75940W | 260 new 
KT998966 | 
LK Saline River, AR 34.58710N, 92.60480W 270 new 
KT988968 
KT988953 | Big Creek, MO 37.29375N, 90.62775W new 
KT988927 _| Big River, MO 37.81483N, 90.77049W M3 new 
KT988941 
Brenton Creek, MO 37.93847N, 90.79239W M4 new 
KT988952 
KT988921 
KT988938 | Meramac River, MO 37.57011N, 91.30293W M5 new 
KT988942 
KT988918 
KT988919 _| Trek Creek, MO 37.95693N, 91.89561W M6 new 
KT988925 
KT988949 | 
Little Piney River, MO 37.91063N, 91.90381W M7 new 
KT988960 
KT988924 | 
LK Mill Creek, MO 37.87296N, 91.92921W M8 new 
KT988928 
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Taxon Accession | Locality Identifier | Reference 
KT988947 
Clifty Creek, MO 38.04173N, 91.96089W M10 new 
KT988948 
KT988914 
KT988915 | North ForkWhite River, MO 36.66723N, 92.28129W M14 new 
KT988926 
KT988931 | James River, MO 37.19040N, 93.12660W M16 new 
KT988937 | Indian Creek, MO 36.79320N, 94.24380W M19 new 
KT988930 | Spring River, MO 37.11560N, 93.89420W M20 new 
KT988961 
Clear Creek, MO 37.30831N, 93.50060W M22 new 
E. potosiensis KT988970 
KT988959 ; 
Niangua River, MO 37.51970N, 92.98420W M24 new 
KT988969 
KT988917 ple 
Big Piney River, MO 37.32720N, 92.00210W M27 new 
KT988946 
KT988935 
RTORRORE Current River, MO 37.27990N, 91.40600W M31 new 
KT988912 
KT988913 | Sallisaw Creek, OK 35.57660N, 94.83047W OK new 
KT988920 
E. showalteri U73767 Lydeard et al. 1997 
E. virginica 1 DQ311117 Lee et al. 2006 
E. virginica 2 AF100989 Holzangel and Lydeard 2000 
To fluvialis AF100999 Holzangel and Lydeard 2000 
Juga plicifera AF101004 Holzangel and Lydeard 2000 
Leptoxis ampla 1 U73768 Lydeard et al. 1997 
Leptoxis ampla 2 KF680604 unpublished 
Le. crassa anthonyi AF101001 Holzangel and Lydeard 2000 
Le. dilatata DQ311122 Lee et al. 2006 
Le. foremani KF680592 unpublished 
Le. picta 1 KF680596 unpublished 
Le. picta 2 U73769 Lydeard et al. 1997 
Le. plicata U73770 Lydeard et al. 1997 
Le. praerosa AF101002 Holzangel and Lydeard 2000 
Le. taeniata 1 U73771 Lydeard et al. 1997 
Le. taeniata 2 KF680600 unpublished 
Le. virgata AF101000 Holzangel and Lydeard 2000 
Lithasia armigera AF100998 Holzangel and Lydeard 2000 
Li. duttoniana AF100997 Holzangel and Lydeard 2000 
Li. geniculata fuliginosa | AF 100996 Holzangel and Lydeard 2000 
Li. geniculata geniculata | AF100995 Holzangel and Lydeard 2000 
Pleurocera acuta | AF100994 Holzangel and Lydeard 2000 
P. acuta 2 MF357697 | Mulberry Fork Little Red River, AR 35.74210N, 92.33380W new 
P. acuta 3 MF357698 | Warm Fork Spring, AR 36.49580N, 91.53320W new 
P. annulifera W777 Lydeard et al. 1997 
P. prasinata | U73774 Lydeard et al. 1997 
P. prasinata 2 U73773 Lydeard et al. 1997 
P. pyrenella 1 AF100990 Holzangel and Lydeard 2000 
P. pyrenella 2 DQ311123 Lee et al. 2006 
P. pyrenella 3 KT164352 Whelan and Strong 2016 
P. uncialis hastata AF100993 Holzangel and Lydeard 2000 
P. vestita U73775 Holzangel and Lydeard 2000 
P. walkeri AF100992 Holzangel and Lydeard 2000 


were clear and reliably scored in each individual. Bands 
were identified and sized using GelAnalyzer2010 (http:// 
www.gelanalyzer.com/), then binned into markers with 
5 bp tolerance, For both studies bands were scored as 
present or absent in each individual to produce separate 
data matrices. 


We analyzed our data matrices using GenAIEx 6.052 
(Peakall and Smouse 2006) and the ‘poppr’ 2.2.1 pack- 
age (Kamvar et al. 2014, Kamvar et al. 2015) in R 3.3.1 
(R Core Team 2016). For each, we generated a genotype 
accumulation curve (Kamvar et al. 2015) to determine 
the minimum number of loci necessary to discriminate 
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Figure 2. Map showing collecting localities for both ISSR studies. Inset shows position of Oklahoma (OK) relative to Arkansas 
(ARK1-3) populations. Shell of E. potosiensis from population OK is shown. 


Table 2. Sequences, annealing temperatures (°C), and number of bands produced for ISSR primers used in each study. 


‘small’ study ‘large study’ 
Primer sequence | Annealing temperature | number of bands | Primer sequence | Annealing temperature | number of bands 
5'-(AG),T 50° Pp (AC), C Bs” 47 
5'-AC, 49° 27 
5’-BHB(GA), 49° 57-(CA),RG 53° 26 
B'-RY(CA), 49° ei et 
5’-WB(GACA), 50? 4 


between individuals in a population. We calculated Simp- 
son diversity and evenness (Gritinwald et al. 2003) at 
each locus, and Nei’s gene diversity (Nei 1973) for each 
population. An analysis of molecular variance (Excoffier 
et al. 1992) was performed on the data by populations 
in regions (spring [S] versus creek upstream [C] versus 
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creek downstream [B], or Arkansas versus Oklahoma) 
to partition the overall genetic variation. The ® statistic 
was used to describe ratios of between and among popu- 
lation and region variance to total variance; © values are 
comparable to F statistics calculated for dominant-only 
genotype markers (Excoffier et al. 1992). Significance of 
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@® was determined through standard permutation (10,000 
replicates) of the entire dataset. We also used principal 
coordinates analysis (PCoA) based on Nei’s genetic dis- 
tance to visualize population structure. 


Results 


Broad scale phylogenetic analyses 


Within EF. potosiensis we observed nucleotide diversity 
m = 0.066 and Tajima’s D = 0.894 (p [D > 0.894] = 0.38) 
using aligned sequences prior to processing with Gblocks. 
Our parsimony network (Figure 3) connected the 30 re- 
sulting 16S haplotypes with a maximum of 83 mutations. 
Our final data matrix processed through Gblocks consisted 
of 110 sequences and 317 nucleotide sites. Under maxi- 
mum likelihood we generated a single tree (log likelihood 
= -2454.1915, se = 186.2554) under the HK Y+R3 model 
(Figure 4). The E. potosiensis sequences did not form a sin- 
gle clade but rather resolved in four well-supported clades. 
Constraining the monophyly of E. potosiensis, however, did 
not produce a significantly different topology (p = 0.63). In- 
dividuals from single localities were not necessarily recov- 
ered as sister taxa, and constraining them to be sister pro- 
duced a significantly less likely topology (p < 0.001). Also, 
no grouping of sequences by river drainage was observed 
in the overall tree nor in the four EF. potosiensis clades, and 
constraining the monophyly of each river drainage generat- 
ed significantly less likely topologies (p < 0.001). 


10 samples 


1 sample 
Arkansas 
Meramac 
St. Francis 
Missouri 
Quachita 


Red 


Figure 3. TCS parsimony network of FE. potosiensis 16S haplo- 
types color coded by river drainage. Hash marks represent mu- 
tational steps between haplotypes. 
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Figure 4. Maximum-likelihood phylogeny of EF. potosiensis 
sequences. Individual snails are color-coded by river drainage 
and labelled according to Figure 1 and Table 1. Branches with 
bootstrap support > 80% are labeled. 
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Population genetics 


In the ‘small’ study, we were able to bin and identify 22 
markers. There were 110 unique genotypes identified 
among the 120 total individual snails. A genotype accumu- 
lation curve suggested that the minimum number of mark- 
ers needed to discriminate individuals was 19, so our anal- 
yses were close to the detection limit for those markers. 
Simpson diversity for individual loci ranged from 0.033 to 
0.500, while evenness ranged from 0.383 to 1. Nei’s gene 
diversity by population ranged from 0.248 (B3) to 0.394 
(C3). AMOVA (Table 3) showed that most of the variation 
in the data was within populations (95.86%), followed by 
between the regions (spring and creek sites) (3.91%) and 
within each region (0.23%). Permutation tests suggested 
significant (p < 0.01) genetic structuring between popula- 
tions and between regions (spring or creek), but not be- 
tween populations within the same region (p = 0.28). Pair- 
wise comparisons of populations suggested low genetic 
differentiation between populations (® < 0.133), and PCoA 
showed no clustering by proximity or region (Figure 5). 


coordinate 2 


In the ‘large’ study, we were able to bin and identify 
a total of 100 markers. Each of the 200 snails possessed 
its own unique genotype. A genotype accumulation curve 
suggested that a minimum of 31 polymorphic markers was 
needed to discriminate individuals, so our dataset of 98 
was sufficient for further analysis. Simpson diversity for 
individual loci ranged from 0.303 to 0.500, while evenness 
ranged from 0.376 to 1. Nei’s gene diversity for the pop- 
ulations ranged from 0.254 for ARK1 to 0.328 for ARK3. 
An AMOVA (Table 4) showed that most of the variation 
in the dataset was seen between individuals of a popula- 
tion (66.8% of total). The least variation (12.1%) was seen 
between populations in a region (Arkansas or Oklahoma), 
and the remainder (21.1%) existed between regions. Per- 
mutation tests indicated significant genetic structuring 
within populations, between populations, and between 
regions (® = 0.133 to 0.367); pairwise comparisons at all 
levels were significant at p<0.01. PCoA suggested com- 
plete separation of OK from the ARK populations, and lit- 
tle overlap between the three ARK populations (Figure 6). 


coordinate 1 


Figure 5. Principal coordinates analysis of ISSR genetic distances from the ‘small’ study showing overlap of populations. Popula- 


tions labeled as in Figure 2. 
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coordinate 1 


Figure 6. Principal coordinates analysis of ISSR genetic distances from the ‘large’ study showing separation of populations. Popu- 


lations labeled as in Figure 2. 


Table 3. Results of the ‘small’ study AMOVA. Regions were 
spring (S populations), and creek above (C) and below (B) con- 
fluence with the spring. 


Source d.f. @ % total 
variation | value 

Within populations 108|0.054} 95.86 |0.014 

Between populations within 10 10.007! 0.23 10.274 

regions 

Between regions 1 | 0.047 3.91 0.005 

Discussion 


Freshwater taxa from the same region tend to share evo- 
lutionary history, resulting in replicated patterns of bio- 
geography and speciation (Walker and Avise 1998, Carini 
and Hughes 2006). This is true in the interior highlands, 
where the Ozark and Ouachita drainages represent sep- 
arate areas of endemism (reviewed in Hoagstrom et al. 
2014). Phylogenies of freshwater mussels reflect these 
patterns (Serb and Barnhart 2008, Inoue et al. 2013, 
Chong et al. 2016), but they remained unexplored in 
freshwater snails. Our nuclear ISSR data suggest that 
Elimia potosiensis exhibits genetic population differen- 


Table 4. Results of the ‘large’ study AMOVA. Regions were 


either Arkansas or Oklahoma. 
% total p 
variation | value 


<0.01 


} 2 Joss 12M Oy |) -aQ70s1 


<0.01 


Source 


Within populations 
Between populations within 
regions 

Between regions 


tiation, even at small geographic scales, while our mi- 
tochondrial data show no evidence of phylogeography 
across the species range. 

Our analyses of mitochondrial 16S sequences suggest- 
ed no genetic structuring. Mitochondrial haplotypes did 
not group together in our TCS network, and we recovered 
E. potosiensis in four well-supported clades instead of 
a single clade. However, we could not reject the mono- 
phyly of E. potosiensis based on our data. No statistical 
support existed for genetic structure within populations, 
within river drainages, or within populations or drainages 
in each of the four clades. We observed no evidence of 
heteroplasmy to suggest nuclear analogs of mitochon- 
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drial sequences (“numts”) in our data, nor did we find 
any biogeographic groupings of sequences that would 
suggest cryptic taxa. What we did observe, however, was 
a pattern where E. potosiensis sequences comprised mul- 
tiple well-supported divergent clades. As mentioned pre- 
viously, this is frequently observed in freshwater snails 
including pleurocerids using mitochondrial markers. 
Whelan and Strong (2016) showed that the pattern ex- 
ists across multiple valid species, is not limited to spe- 
cific taxa, and becomes more apparent as additional se- 
quences from each population are included in an analysis. 
The maximum number of individuals we sequenced for 
a single population was five, and two of these individu- 
als exhibited the highest pairwise sequence divergence. 
We predict that additional EF. potosiensis samples would 
only increase the observed sequence diversity while not 
increasing phylogenetic resolution. Whelan and Strong 
(2016) proposed that balancing selection might keep di- 
vergent haplotypes in pleurocerid populations. Our Ta- 
jima’s D of 0.894 suggests balancing selection of 16S 
haplotypes in EF. potosiensis, though population reduction 
and migration can generate similar results (Maruyama 
and Fuerst 1985, Simonsen et al. 1995). 

Four nominal morphospecies comprise the modern no- 
tion of E. potosiensis (Burch and Tottenham 1980): E. po- 
tosiensis (Lea, 1841), a narrow-range species from Mis- 
sourl; E. crandalli (Pilsbry, 1890), endemic to Mammoth 
Springs, Arkansas; E. ozarkensis (Call, 1886), found in 
springs in north Arkansas; and EF. plebejus (Gould, 1850), 
the most widespread and the shell form currently associ- 
ated with EF. potosiensis. While E. potosiensis does pos- 
sess shell variation that may have an environmental com- 
ponent (Minton et al. 2011), no additional data supports 
recognition of any separate morphological entities. Jones 
and Branson performed detailed examinations of internal 
anatomy and radula structure of all four nominal E. po- 
tosiensis taxa, noting that, without their shells, internal 
structures of the morphotypes “...can not be separated” 
(Jones and Branson 1964: 60), and that the shell varia- 
tions seen among the forms could be found throughout 
the species’ range. They thus treated EF. potosiensis as 
single species exhibiting shell plasticity, where morpho- 
logical “...variation is the rule rather than the exception” 
(Jones and Branson 1964: 60). This anatomical data, 
combined with our inability to reject the monophyly of E. 
potosiensis, lead us to support the current notion (Burch 
and Tottenham 1980, Johnson et al. 2013) that E. poto- 
siensis is a Single widespread species. 

While we saw no broad scale genetic patterns in E. po- 
tosiensis, our two ISSR studies suggested population-lev- 
el genetic structuring. Our ‘small’ study of ten individuals 
from each of 12 populations represented a pilot study as 
well as our first effort with ISSRs. In it, we were able to 
characterize 22 bands visualized on agarose gels using six 
primers. The AMOVA suggested that most of the genet- 
ic variation was within each population, but that a small 
(Hartl and Clark 1997) yet significant amount of differ- 
ence existed between the spring and two creek regions. 
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We predict that the spring 1s isolated from the creek by 
the one meter of elevation, as the spring run flows over 
a rocky ledge and into the creek below. This normally 
serves as an isolating mechanism between populations in 
the two channels. Snails within each channel can migrate 
upstream and downstream freely, and during heavy rains 
and flooding conditions, the creek can crest above the 
spring. This could potentially carry individuals from the 
spring into the creek and vice-versa, and may explain the 
low pairwise ® values observed between populations and 
lack of PCoA clustering (Wilmer et al. 2011). Our ‘large’ 
study showed complete differentiation of the four popu- 
lations examined with high levels of genetic differentia- 
tion at all population and region levels (Hartl and Clark 
1997). We used large numbers of individuals and mark- 
ers, both of which increase the resolution and accuracy of 
population genetic studies (Ruzzante 1998, Kalinowski 
2005), and utilized polyacrylamide gels to better visual- 
ize bands. Our AMOVA and PCoA suggested that each 
population is genetically distinct from the others. 

Two issues we identified in our ‘small’ study were the 
low number of bands resolved by our primers and the 
small sample sizes from each population. Our band total 
for six primers was comparable to the number generated 
by a single primer in other studies, and is likely indicative 
of two factors. First, choice of our primer sequences was 
probably not ideal, employing dinucleotide repeats and/ 
or degenerate 5’-bases ahead of the repeated portion. Re- 
sults from the literature suggest that trinucleotide repeats 
and 3’-degenerate anchors after the repeats increase res- 
olution and repeatability (Godwin et al. 1997). Second, 
we visualized our bands on agarose instead of polyacryl- 
amide. Agarose has far less resolving power, and multiple 
bands in a small range of sizes might not have been iden- 
tifiable on our gels (Stift et al. 2003). In regards to sam- 
ple size, many studies recommend thirty individuals per 
population (Pruett and Winker 2008, Hale et al. 2012). 
Our AMOVA suggested treating the regional groupings 
as populations, where we were able to resolve significant 
structuring on a small spatial scale given the now larger 
sample sizes. The differences in our ‘small’ and ‘large’ 
studies highlighted the importance of having the appro- 
priate number of markers and individuals. 

Our study also highlighted the need for novel and use- 
ful genetic markers for pleurocerid conservation. Near- 
ly 80% of all pleurocerids are imperiled (Johnson et al. 
2013), and most species designations still rely on shell 
morphology. The genetic species and population structure 
of most pleurocerid taxa remains unknown. Our study 
provides the latest example where gene sequences are not 
useful in delineating neither population nor species-lev- 
el boundaries. Results from analyses of mitochondrial 
gene sequences are mixed, generating far less resolution 
(e.g. Minton and Savarese 2005, Minton 2013) than un- 
certainty (reviewed in Whelan and Strong 2016). Single 
mitochondrial gene sequences, even from the same popu- 
lation, can differ by over 20%, which can lead to hypothe- 
ses of species non-monophyly and aberrant biogeography 
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(Dillon and Robinson 2009). Nuclear gene sequences, on 
the other hand, are too highly conserved in pleurocerids 
to be useful (Whelan and Strong 2016). Other forms of 
nuclear data, however, such as allozymes and now ISSRs, 
have demonstrated utility in population genetic studies. 
Allozymes have been used to show genetic connectivi- 
ty in various pleurocerid genera (Chambers 1980, Dillon 
1984, Stein and Stansbery 1984, Dillon 2014), and herein 
our data suggests ISSRs show promise across large and 
small geographic levels. No published applications of 
more modern population genetic techniques (microsatel- 
lites, SNPs, etc.) are available, but we hope that they too 
will be useful in understanding pleurocerids. The contin- 
ued development of consistently useful tools is needed 
to help properly understand genetic population and spe- 
cies structure in pleurocerids, especially in a conservation 
contextS before more taxa are lost. 
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